Strategy for implementing stabilizer-based codes on solid-state qubits 
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We present a method for implementing stabilizer-based codes with encoding schemes of the opera- 
tor quantum error correction paradigm, e.g., the "standard" five-qubit and CSS codes, on sohd-state 
qubits with Ising or XY-type interactions. Using pulse sequences, we show how to induce the effec- 
tive dynamics of the stabilizer Hamiltonian, the sum of an appropriate set of stabilizer operators for 
a given code. Within this approach, the encoded states (ground states of the stabilizer Ifamiltonian) 
can be prepared without measurements and preserved against both the time evolution governed by 
the original qubit Hamiltonian, and energy-nonconserving errors caused by the environment. 

PACS numbers: 03.67.Lx,03.67.Pp,03.67.Ac 



A variety of quantum error-correcting codes (QECC) 
have been widely investigated aiming at a robust comput- 
ing system similar to the classical digital computer [l]- 
In particular, codes based on the stabilizer for- 
malism constitute an important class of QECCs. This 
formalism has proven useful not only for the standard 
codes [ll 2|[, but also for the subsystem code ^tM^ topo- 
logical^ ^fl 15 1, and Majorana codes @. On the ex- 
perimental side, Knill et al. demonstrated its usefulness 
in the NMR domain [Uliii. Stabihzer-based QECC in 
systems with always-on coupling have recently attracted 
a great deal of interest [l^ Il3| • 

Stabilizer operators Gj {j = 1,...,/) are mutually 
commuting operators given by products of multiple Pauli 
matrices X^, Yi, and (i = 1, . . . , n) [2]. Conventionally, 
logical qubit states are encoded through measurements 
into a joint, 2'-dimensional, eigenspace 'Hs of these op- 
erators. For I stabilizer operators and n physical qubits, 
a maximum number oi k = n — I logical qubits can be 
encoded into H5, while k < n — I \n case of subsys- 
tem encoding. Since the ground states of the stabilizer 
Hamiltonian iTstab '■= — ^3 ^^^^ joint eigenstates 

of all stabilizer operators, its ground state manifold can 
play the role oiHs- 

It is important to note that stabilizer operators of 
many error-correction codes, e.g., the surface code Q 
or color code , are given by products of more than two 
Pauli matrices. Therefore the corresponding stabilizer 
Hamiltonians cannot be directly implemented in natural 
solid-state qubit systems, where the interactions between 
qubits are of two-body type [H, [l^ . 

In this letter, we demonstrate how to prepare ground 
states of Hstah as encoded states and preserve them by 
inducing the effective dynamics of this Hamiltonian using 
sequences of pulses in the form of singie-qubit rotations. 
Being based on single-qubit rotations only, our method 
works for an always-on physical (qubit) Hamiltonian, i.e., 
it does not require switching on and off any of its parts 
(single-qubit or interaction). 

The distinguishing feature of this method is that it al- 



lows the preparation of encoded states without measure- 
ments, thus avoiding measurement-induced decoherence. 
The method can be used not only for standard codes (i.e., 
five-qubit and CSS codes) but also for the extended class 
of codes with encoding schemes within the general op- 
erator quantum error correction framework [j, Q . Even 
in the presence of inevitable pulse (rotation angle) errors 
the ground-state fidelity scales favorably with the system 
size. 

Our scheme allows the implementation of stable solid- 
state quantum memories. ThiSj in turn, facilitates the 
realization of quantum gates |l7| within the limitations 
imposed by the size and coherence time of a system. 
Although the scheme requires a rather large number 
of pulses (rotations), its feasibility can be anticipated 
based on the recent progress in qubit-manipulation tech- 
niques [id ]. 

Generation of stabilizer operators. — We consider the 
Hamiltonian H ~ Hq + Hxy, where Hq — ^^Hoi — 
J2i{^iXi + £iZi) is a single-qubit part and Hxv — 
H'^+^ = J,- ,+i(X,X,+i +y,r,+i) the (two-body) 
Xy-interaction part. Note that instead of the XY- 
Hamiltonian we could also use the Ising Hamiltonian 

-falsing = X/z Ji,i+lZiZi+l- 

Higher-order products of Pauli matrices can be gener- 
ated using the following transformations 2l|: 



^xY Yie^ 



cgXi — seZiYij^i , 

CgY^ + SgZiXi+i , 



+ cese{X,Y,+i^Y,X,+i), (1) 

where eg = cos(26') and se = sm{29). For 6 = Jt = tt/A, 
these transformations increase the order of the Pauli- 
matrix terms as Xi —ZiYi-^.i and Yi — ZiXi^i. Sim- 
ilarly, one obtains Zi — )> Zi^i. Analogous relations hold 
for the Ising interaction. 

Manipulation of the Hamiltonian. — The key step in 
obtaining the local stabilizer operator is extracting a 
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process : e '^opHop jc^giTopHop 


time for generation 


Gi 
Ga 
G3 
G4 


X Z Z X I 
I X Z Z X 
X I X z z 
Z X I X z 


^2X2 

02X2 

^2X2 


J/fy ^ (^/2)i ^ /^iV + ^^iV 

Gi iffy ^ (7r/2)i(7r/2)i ^ Hfy 

Gi ^ H*!y ^ (7r/2)i(7r/2)g ^ i/fy ^ (^/2)?(7r/2)| 


Tini + 24rrot + 4rop 
Tini + 24rrot + 4rop 
Tini + 43rrot + STop 
Tini + 45rrot + 8rop 



TABLE I: Stabilizer operators for the five-qubit code 1] as realized with XF-type interactions. Each appearance of H^xy^ 
indicates an application of the transformation H g-^^opifop j^g»Top-tfop^ while each instance of (7r/2)f denotes a 7r/2 rotation 
about the a;-axis. The rightmost column shows the time required to generate each stabilizer operator; Trot is the time needed 
to perform a single-qubit rotation. 



single-qubit part or a pure two-body interaction part 
from the original system Hamiltonian. This process is 
carried out using the Baker-Campbell-Hausdorff (BCH) 
formula [22*1 . For simplicity, we explain this procedure for 
ffi = in iJo, where only rotations about the z axis will 
be needed, and set Vli — Vl. In the general case, the pro- 
cedure requires a slightly more complex pulse sequence. 

A part Ha can be extracted from Hq by applying a 
single appropriate tt pulse, if that pulse transforms to 
Ha — Hb, where Hh = Hq — Ha consists of the unwanted 
terms. For 2n alternating periods of propagation with 
A = iT{Ha+Hb) and B = iT{Ha-Hb), the BCH formula 
yields 



(e^e^)" « exp(i2nTiJ, + nT'[Ha,Hb]) 



(2) 



where the duration of the pulse sequence is 2nT. Thus, 
as long as T\\Hb\\ < 1, where \\A\\ = \It:{A'^ A) / d]'^/'^ is 
the standard operator norm in a Hilbert space of dimen- 
sion d, we can neglect the second term. As the number 
n of repetitions increases, this approximation becomes 
progressively better. 

In order to extract a single-qubit (local) part of the 
system Hamiltonian, relation ^ has to be applied twice, 
leading to (case n = 1) 

e^e^e^'e^' « exp(2/i, + [hb, ha]) exp(2/i; - [h^ h'a]) 

« exp(2(;^, + h'a) + [hb, ha] - [hi h'a] + ^[ha, h'a])i?,) 



where h^'\ := irH^''^ 



A' = K+ /i' and B' = K- h',. 



a/b ■ a/b' '^a''% "^^^^ ^ "'a 

Consider extracting a single-qubit Hamiltonian X2 for 
a one-dimensional five-qubit array. With hi = irflXi 
(single-qubit Hamiltonian with = 0) and hij = itH'^yi 
we set 



ha = 

hb = 

h'a = 

K = 



h2 + h34 + ft.45, 

hi + /1.3 + /l4 + /l5 + ^112 + ^23, 
^2 — ^34 — ^45 1 

hi - h^ + hi - h^ + hi2 - ft.23- 



(4) 



The sequence of operators describing the time evolu- 
tion on the left-hand side of Eq. ([3]) is obtained in the 
following manner: by applying a 7r-pulse to qubits 1, 
3, 4, and 5 one transforms A into B and B' into A', 



while B is transformed into B' by a 7r-pulse applied 

'ini — JJ^v^^ 



to qubits 3 and 5. This leads to 



where H-^-^^ = ^1X2 is the desired initial Hamiltonian and 

H^^l = [r!r/2] (732^3^2 - JsiYsZi - J^^Y^Z^) is an un- 
wanted perturbation term. This term scales like Jy-r <C 1 
and hence can be reduced by shortening the duration of 
the pulse sequence. 

Similarly, the operator e~^^^^Y is obtained by extract- 
ing H^Y from the system Hamiltonian using 

ha = h23 + hi + hi, 

hb = hi2 + /l34 -I- ft.45 + h2 + /l3 + /l5, 

h'a^ h23 - hi~ hi. 



mi 



h'^- 



-hi2 



'34 



'45 



h. 



(5) 



The perturbation terms can be neglected for J/Q. <^ 1. 

Time- evolution of the system — To illustrate the time 
evolution of the system during the generation process for 
the stabilizer operator Gj we will use the schematic nota- 
tion p(0) pit). Here p{t) = exp{-iHt)p{0) exp{iHt) 
is the density matrix for a time-independent Hamilto- 
nian H, or for an effective H in the sense of average- 
Hamiltonian theory (2^ . The manipulation of the orig- 
inal Hamiltonian iJini is described by [2l[ p(0) 
r,niH^n, jjp^op ^^^^^ where Top = 7r/(4J) and Hop is 
either the pairwise interaction H^xy ^ generator of a 
single-qubit rotation, see Table I. During this pulse se- 
quence the stabilizer operator Gj describes the dynamics 
of the system in the sense of average-Hamiltonian the- 
ory 2l|. Here Tjni has to be chosen such that the entire 
process can be carried out in a time interval sufficiently 
shorter than the coherence time. 

In the following we apply our scheme to the standard 
codes (five-qubit and CSS codes) , as well as the surface 
code [3|]. 

Realization of the five-qubit and CSS codes. — The gen- 
eration processes of the four stabilizer operators Gj (j = 
1 , . . . , 4) of the five-qubit code [H are shown in Table 
I. For example, starting from the initial Hamiltonian 



H, 



^2X2, the stabilizer operator Gi of the five-qubit 
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code is realized through the sequence 

g-»ropHjVx2e*''°P^^'^ ^ -Z2Y3 
_g-»(^/4)X.^^^^g»(^/4)X2 ^ Y2Y3 (6) 

g-»r„p[Hi-V+-ffiV]y2r3e'''°p[^i'^+^i''^l ^ . 

The minimal time required for this process is Ti^i + 
24Ti-ot + 4rop. The effective dynamics of i?stab = 
— X^iLi is induced by subsequent generation of the 
four stabilizer operators. 

We would now like to address the feasibility of this 
scheme in a typical superconducting qubit system. For 
two superconducting qubits in a circuit-QED setup the 
resulting effective intcr-qubit interaction is also of XY 
type [11 m. For instance, for g/A = 0.1, g/(27r) = 
200 MHz, A/(27r) = 2 GHz, where g is the Jaynes- 
Cummings coupling constant and A the detuning be- 
tween the resonator frequency and the qubit split ting , we 
have J/(27r) — 20 MHz. Assuming Trot ~ 1 ns [25|, we 
obtain a minimal total time of r^ojc = 24rop + 136ri.ot ~ 
300 ns, which is significantly shorter than T2 ^ 20 /iS 
reported in [26{ . 

Table H shows how to generate the CSS code. 
The stabilizer operators G4, G5, Gg are obtained by 
e-^i:,Y^/i[Gi + G2 + G3)e''^^^^/^. Thus, the minimal 
total time is r™gg = 44Top + 246Ti.ot ~ 600 ns, i.e., again 
much shorter than the T2 given in Ref. j2Q] . 

Realization of the surface code. — Four types of stabi- 
lizer operators are needed to realize the Hamiltonian of 
this code. Qubits are placed at the edges of the square 
lattice, see Fig.[T] Stabilizer operators Hg — ^jestaiis)^j 
are assigned to each vertex s, and Hp — Tlj ehoundary (p)Zj 
to each face p. Using the relations Xi —ZiYi^i 
and Yi — > ZiXi^i, we can form products of nearest- 
neighbor operators such as Yi — ^ Z1X2 — >■ —Z1Z2Y3 ^■ 
In generating adjacent stabilizer operators, 
care should be taken to avoid mixing them. This can 
be achieved by decompositions like Hs — Hs^ + Hs^ [see 
Fig. [Ha) and (b)] and Hp = Hp, + Hp^ [Fig. [UJc) and 
(d)]. The four elements can then be combined into the 
total surface-code Hamiltonian. 

Preparation of encoded states — Our approach also al- 
lows us to prepare encoded states (or codewords) of gen- 
eral stabilizer-based codes without performing measure- 
ments on the system. We show this in detail for the stan- 
dard codes, which encode k logical qubits into a subspace 
of dimension 2^^. However, this procedure also works for 
subsystem encoding provided suitable stabilizer opera- 
tors are added. For any given code, only those Gj with 
1 < J < m and m < n — k that contain X oi Y operators 
are needed for the preparation: 

\ci...ck) = (l + Gi)---(l + G„,)Xi^^---X,^'=|0...0) 

k 71% 




FIG. 1: Four types of surface-code generation on the 2x3 
lattice in Ref. starting from a single-qubit Hamiltonian 
that includes Y operators. All operators in (a), (b) are trans- 
formed to X and those in (c), (d) to Z. Subsequently, the four 
types of stabilizer operators are combined into the topological 
Hamiltonian in the manner discussed in the text. 

where c,; =0,1 and operators Xi act in the logical state 
space {|0)i, |T),;}. Here, Gj denotes a modified stabilizer 
operator obtained from Gj by replacing one of the X 
operators — acting on, say, qubit a — by a Y operator, or 
vice versa. This is done in order to match the effect of 
exp[i(7r/4)Gj] with the action of the projector (1 -f Gj) 
when qubit a is in state |0}. 

We illustrate the encoding procedure on the example of 
a three-qubit code whose stabilizer operators are ^1X2 
and X2X3. This is realized in a three-qubit system with 
Ising interactions. The stabilizer Hamiltonian J{XiX2-\- 
X2X3) is obtained by removing the single-qubit part Hq 
of the original Hamiltonian using 7r-pulses. Its ground 
states are the degenerate states |0} = |000) + |110) + 
|011)-|-|101) and |1} = |001) + |010) + |100) + |111). They 
are obtained by |0)l = \{l + A:iA:2)(1 + A:2A:3)|000) = 
exp[i(7r/4)YiX2]exp[i(7r/4)A:2r3]|000). 

For the five-qubit code, we have Gi = 1^1^2^3X4, 
G2 = X2Z3ZiY^^ G3 — X1Y3Z4Z5, Gi ~ Z1Y2X4ZZ, 
and the multiplication in Eq. (jT]) is carried out in the fol- 
lowing order: exp[z(7r/4)G2] exp[i(7r/4)G4] exp[i(7r/4)G3] 
exp[i(7r/4)Gi]_. In the case of the CSS code, Gi = 
XiX2X3Yi, G2 = XiX2X5Ye, G3 = X^X^X^Yr. Note 
that here only three out of six stabilizer operators are 
needed for the preparation of an encoded state. 

Robustness against pulse errors. — Since the codeword 
states are encoded in the twofold-degenerate ground- 
state manifold |0) and |1) of i/stab, the robustness of 
this method is limited by the rate of leakage out of this 
manifold. In principle, precise estimates of the leakage 
due to the thermal environment could be obtained by 
studying the stability of the ground state to various per- 
turbations as in Ref. However, already single-qubit 
errors — in most cases the prevalent kind of errors created 
by a thermal bath — are suppressed by a factor of e~^^/^ 
or smaller for the five-qubit and CSS code [11]. Hence, 
at low temperatures, unavoidable pulse imperfections are 
likely to be the predominant cause of leakage. 

To estimate this effect, we consider pulse errors that 
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process : e '^op^op jy-giTopHop 



iV ^ (^/2)i(^/2) 

rr23 i rr45 
0-23 , rr45 



time for generation 



Gi 
G2 

G3 



X X X X 7 7 7 
X X 7 7 X X 7 

X 7 X 7 X 7 X 



^2X2 
-^3X3 

-t^3X3 



77-34 

(7r/2)I(7r/2)4-(7r/2) 

77iV ^ (^/2)f ^ 77iV + 
(7r/2)^(7r/2)^(7r/2)^(7r/2)g 



rrl2 
tLxY 



'xV^(vr/2)| 

rr23 , rr45 > rrl2 , 0-56 
nXY + J^XY ~^ ilXY + J^XY ~ 
ri-12 , rr34 , u-56 
J^XY + ^XY + ^XY 



J^XY 



Tr67 
• ^XY 



Tini + 26rrot + 4rop 
Tini + 45rrot + STop 

Tini + 51rrot + IOTop 



TABLE II: Stabilizer operators for the CSS code 01 as realized with XY-type interactions. The operators G4, G5, and Ge are 
obtained by replacing X with Z in Gi, G2, and G3, respectively. The rightmost column shows the time required to generate 
each stabilizer operator; Trot is the time needed to perform a single-qubit rotation. 



can be modeled by randomly distributed, unbiased, and 
uncorrelated deviations 66 vi^ith ag = from the 

ideal angle of tt/2. The leakage can then be estimated 
by looking at the average of the ground state fidelity 
F{t) = |(0|J7p(t)|0)p, where Up{t) is the time evolution 
operator with imperfect pulses. This average is approxi- 
mately given by {F{t)} « 1 - Npa'^t/{8T), where A'^p is 
the number of pulses in the sequence to generate TJstab, 
and T its duration. The number iVp of pulses is given by 
the number of rotations needed to generate all stabiliz- 
ers of a given code (for the five-qubit and CSS code, see 
Table Hand Table III respectively). 

In conclusion, we have demonstrated measurement- 
free preparation of encoded states in the stabilizer-based 
codes described by the operator quantum error correc- 
tion paradigm. The scheme is based on pulse sequences 
applied to solid-state qubit Hamiltonians with two-body 
interactions of XY- or Ising-type. We have estimated the 
intrinsic robustness of our scheme against pulse imperfec- 
tions. Once implemented experimentally, our scheme will 
pave the way for robust quantum information processing. 

We would like to thank Fabio Pedrocchi and Adrian 
Hutter for discussions. This work was financially sup- 
ported the Swiss SNF and the NCCR Quantum Science 
and Technology. 
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